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The rapidly increasing integrated luminosity collected by the CDF and DO detectors at the Tevatron Collider has 
qq ■ resulted in a wealth of new results from searches for Higgs bosons in the extensions of the Standard Model. Tighter 

f***^ ' limits are set on the parameters governing the Higgs sector in these models. 

o" 

1. EXTENSIONS OF THE STANDARD MODEL HIGGS SECTOR 

£f-\ | It is well known that the Standard Model (SM) of particle physics, successful as it is in describing the strong, 

electromagnetic, and weak interactions, cannot be complete. It does not allow for the incorporation of gravitation; 
^> . and the mass of its scalar Higgs boson receives quadratic radiative corrections, making it unnatural for it to be as 
1* \ low as the order of the electroweak symmetry breaking scale, several hundred GeV. 

Qh The most popular extension of the SM designed to address these issues is Supersymmetry or SUSY. For each SM 

particle, it introduces a partner with spin differing by half a unit but with otherwise the same properties. This allows 

l — '■ for a cancellation of radiative corrections due to loops of SM particles and their superpartners. But Supersymmetry 

,__! j must be broken, for otherwise the superpartners' masses would equal those of the SM particles, and SUSY would 

have long been discovered. The precise breaking mechanism depends on the interactions at energy scales inaccessible 

2j^ ■ by the present-day colliders, and is not known. 

A practical approach towards a description of SUSY phenomena at the energy scales accessible at experiments 

being conducted or in preparation is to consider the general effective theory obtained upon SUSY breaking: this 

is the Minimal Supersymmetric Standard Model (MSSM). Besides SUSY conserving interactions, it features soft 

^ri ! SUSY breaking interactions that still maintain the cancellation of quadratic corrections to scalar masses, one of the 

reasons to introduce SUSY in the first place. The MSSM constructed in this fashion does have one severe drawback, 
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however: it features 105 new and a priori arbitrary parameters. Even if many of these would give rise to unobserved 
features such as lepton flavour violation, and are hence constrained by measurements, the exploration of the full 
MSSM parameter space is a daunting task. 

In contrast, the same Higgs sector that is at the root of the problems with the SM is described (at tree level) by 
a mere two parameters. SUSY requires the existence of two Higgs doublets (of hypercharge Y = \ and Y = —1). 
Upon electroweak symmetry breaking, this results in five physical Higgs bosons: two charged ones, H ± 1 two CP-even 
neutral ones, h and H (by convention m^ < mjj), and one CP-odd one, A. At tree level, all masses can be expressed 
as a function of only two parameters, usually taken to be ttia and the ratio of the two Higgs doublets' vacuum 
expectation values tan/3 = Ui/va- 

Also the Higgs boson couplings to SM particles can be expressed as a function of tha and tan (3. Table I shows 
the neutral Higgs boson couplings (here, a is not a free parameter but itself depends on tan/3 and m^); the H Uidj 
coupling is proportional to V Ui dj (m Ui cot/3(l — 75) + m^ tan/3(l + 75)). 

A remarkable feature is that (again at tree level) irih < rnz- This sole fact would be enough to exclude the MSSM, 
were it not for the fact that radiative corrections relax these constraints. In principle, these radiative corrections 
involve the whole MSSM, but the number of parameters contributing substantially is still small. As a general matter, 
a constraint m/, < 135GeV remains. 

The radiative corrections affect the Higgs bosons' couplings as well as their masses. Rather than considering the 



Table I: MSSM Neutral Higgs boson couplings to SM particles, relative to the respective SM Higgs boson couplings. 
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whole MSSM parameter space, several representative scenarios are typically considered that evade the bounds set 
by searches for Higgs bosons conducted at the LEP e + e~ Collider. Specifically, the scenarios considered here are 
m™ x scenario (with MSSM parameters tuned to obtain the highest possible mh) and the no-mixing scenario (with 
parameters tuned to switch off the mixing between the t\ and ti superpartners of the top quark) [1]. Within these 
scenarios, the Higgsino mass parameter fi is in principle fixed; however, given the sensitivity to this parameter, both 
/i = +200GeV and /i = — 200GcV are considered in some of the analyses below. 

2. NEUTRAL HIGGS BOSONS 
2.1. Generalities 

Previous searches at LEP for neutral Higgs bosons have excluded a large region of the (m a, tan (3) parameter space 
where the lightest Higgs boson becomes sufficiently light to be accessible at LEP. In particular, the low tan/3 region 
has been excluded up to large values of tua- 

The searches conducted at the Tevatron Collider are largely complementary to those earlier searches. At large 
values of tan/3, the enhanced couplings to b quarks and r leptons lead to significantly increased production cross 
sections. This fact, in combination with the particular signatures of 6-quark jets and r leptons, is exploited below. 

As a consequence of large tan/3, the masses of the A and H bosons (or A and h, for low values of tua) become 
almost degenerate (their mass difference becomes less than the corresponding detector resolution). This allows for 
a significant simplification in the analyses: rather than attempting to identify the two nearby resonances, they are 
simply lumped together, and denoted as <fi below. By comparison, the contribution from the "third" neutral Higgs 
boson becomes negligible. 
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The b quark being the heaviest SM isospin — i fermion, the <j> — > bb branching fraction is the largest one for high 
values of tan/3, becoming as high as ~ 0.9. As a consequence, it constitutes one of the dominant search channels. 
Unfortunately, this decay cannot be identified inclusively, due to an overwhelming QCD bb production background. 
Instead, one looks for the final state in which the cj) is radiated off an initial state 6 quark 1 . A significant (~ tan 2 (3) 
cross section enhancement occurs also in this final state; the identification of a third b jet, typically of lower px, 
suppresses greatly the QCD background. 

The CDF Collaboration has searched for the b(f)(— ► bb) channel in a data sample corresponding to an integrated 
luminosity of 1.9 fb _1 , requiring all jets in a three-jet sample to be 6-tagged by a displaced vertex algorithm [2]. 

Even this triple-tagged sample is dominated by QCD backgrounds. Each jet may be a genuine 6-quark jet, or a 
fake (charm or light-flavour) jet. The two leading jets, hypothesized to result from the <j> decay, are treated on an 
equal footing. Therefore, considering only those final states containing at least two genuine 6-quark jets leads to five 
background categories: bbb, bbc, bbq, beb, and bqb (here, q denotes a light-flavour jet). 



1 Here and in the following, charge conjugated final states are implied as well. 



The distribution of the invariant mass m\i of the two leading jets for the signal is taken from simulation. For the 
backgrounds no adequate simulation exists, and templates for each of the above categories are obtained starting from 
a data sample containing at least two b tags. For the first three categories, this starts from a bbj sample to which 
parametrized efficiencies for tagging the third jet (6, c, or q) are applied; for the bqb category, a similar procedure is 
applied to events with the third and one of the two leading jets already tagged. An additional correction, derived 
from simulation, is applied to account for the different processes contributing to final states with more than two b 
jets. Finally, the bcb distribution is assumed to look like the bbb one, but with only a subset of processes contributing; 
another simulation-derived correction is applied to account for this. 

The normalization of the distributions obtained in this fashion are adjusted to fit the data distribution. However, as 
can be observed from Fig. 1, some of the background templates cannot be distinguished. Therefore, a second variable 
involving the vertex masses (the invariant mass of all tracks, assuming they are n , attached to the displaced vertices) 
is constructed. This is done starting from the same data samples above, and uses the same information in addition 
to known vertex mass distributions for different flavours. 
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Figure 1: Left: invariant mass distributions for various background sources in the CDF 
{mA, tan/3) plane in the m£ lax scenario. 



bb analysis. Right: limits in the 



For each hypothesized value ofm^, two template fits are then carried out jointly to the two distributions, one with 
and one without the assumed signal. No significant improvement in the fit likelihood is observed upon inclusion of 
a signal component; limits on a possible signal component are set in a modified frequentist framework based on [3]. 

A complication is that for sufficiently high values of tan/3, the neutral Higgs boson decay width exceeds the experi- 
mental invariant mass resolution. The resulting broadening of the signal "myi distribution reduces the discrimination 
power between signal and backgrounds, and needs to be taken into account explicitly. 

A similar analysis has been carried out by the DO Collaboration on a 1 fb _1 dataset [4]. An artificial neural network 
6-tagging algorithm exploiting lifetime information was employed. Rather than using a second discriminating variable 
to determine the various background components, multiple fe-tagging criteria (with known efficiencies) are used to 
determine the composition of the double-tagged sample. 

A second discriminating variable, based on kinematic information, is subsequently used nevertheless, to suppress 
further the QCD backgrounds. Two discriminants D are constructed, one combining low-mass (90-130 GeV) and one 
combining high-mass (130-220 GeV) simulated Higgs signal samples. The background predictions are obtained, as 
a function of both the invariant mass and the discriminant, by applying a correction derived from simulation to the 
double-tagged sample. Fig. 2 compares the low-mass discriminant distribution integrated over all invariant masses, 
with its prediction; a good agreement is observed for the background dominated low- D region. The final background 
predictions are obtained by selecting the high-Z? region. The whole analysis is applied jointly to exclusive three-, 
four-, and five-jet data samples. Again no significant excess of data over predictions is observed, and a modified 
frequentist method based on [3] is used to set limits on a possible signal. The results are also shown for both the 
no-mixing and m™ ax scenarios in Fig. 2. A slight excess, not enough to be qualified as a signal, is observed at masses 



between 150 and 200 GeV, leading to a reduced exclusion region compared to the predictions. 
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Figure 2: Left: low-mass kinematic discriminant distribution in the DO <f> — > bb analysis. Right: limits in the (itia, tan/3) plane 
in different scenarios. 
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For high tan/3, the branching fraction for the decay <f> — * t + t~ is approximately 0.1. The drawback of this lower 
branching fraction is offset, however, by the absence of irreducible QCD background. The relative cleanliness of this 
decay channel allows for an inclusive analysis, where the <f> can be produced singly, through gluon fusion gg — ► <p. 

A new DO analysis has been performed on er^, /iTh, and e/i final states (here, only the r decay products are 
indicated and neutrinos omitted; th indicates the r hadronic decay mode) in a 1 fb _1 dataset [5]. Hadronic r decays 
are classified according to their charged track and EM cluster multiplicities, and are selected efficiently using artificial 
neural networks. 

The signal is characterized by two high-px r leptons. To reject W(— ► ^)+jets events with a jet misidentified as 
a Th, the transverse mass reconstructed from lepton and the missing transverse momentum vector (interpreted as 
a single neutrino) must be below 40 GeV in the /iT/j channel and 50 GeV in the er^ channel. Backgrounds from tt 
events are suppressed by vetoeing events with high jet activity. 

The irreducible background from Z — ► t + t~ remains and is estimated using simulation. Signal events are 
recognized as an excess over background predictions in the distribution of the visible mass in v ^ a , computed 
as the invariant mass of the visible r decay products and the missing momentum four-vector approximated as 



i>=[^l$l + $l,$ x ,$, : 
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The visible mass distribution is shown in Fig. 3 for the three channels considered. The expected signal due to a 
160 GeV Higgs boson is also indicated. No significant excess is observed, and limits are again set using the modified 
frequentist framework. It should be pointed out that here, due to the worse experimental resolution caused by the 
escaping neutrinos, the effects of the finite (f> decay width are very modest. 

A very similar CDF analysis has been carried out on a 1.8 fb^ 1 dataset [6]. An earlier analysis, performed on 1 
fb _1 , yielded a small excess that could be attributed to a Higgs boson of mass m^ sw 140GeV. With the increased 
statistics considered here, a good overall agreement between data and background predictions is obtained. This 
analysis now sets the tightest limits in the (tua, tan/3) plane, as displayed in Fig. 4. 



2.4. Fermiophobic Higgs Bosons 

Besides the reasonably "standard" no-mixing and m™ ax scenarios, alternative scenarios are possible where the 
Higgs boson couplings to fermions are small; an example is the "small a e ff" scenario [1], in which radiative corrections 
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Figure 3: Top: visible mass distribution in the three channels considered in the DO 
(niA, tan/3) plane obtained in the m™ ax (left) and no-mixing scenarios. 
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conspire to yield small effective couplings to down- type quarks and charged leptons. 

The DO Collaboration has considered a fermiophobic Higgs model, in which all couplings to fcrmions are suppressed 
but the couplings to gauge bosons are unmodified. In a 2.3 fb -1 dataset, a search was performed for the decay 
h — >• 77 [7]. This would be the main decay mode for relatively light Higgs bosons, rrih < lOOGeV. 

Backgrounds to this process originate from QCD direct 77 events, as well as 7+jet and di-jet events in which one 
or both jets, respectively, are mis-identified as photons. The mis-identification contribution is suppressed using a 
neural network employing variables based on calorimeter shower shapes; the network performance has been measured 
in Z — > e + e~ events and samples of photon-like QCD jets. Starting from a loose preselection, the identification of 
events where neither, one, or two of the photon candidates satisfy a tight neural network cut allows to determine the 
contributions from the individual backgrounds. The irreducible 77 background estimated in this way is consistent 
with the result obtained from a NLO calculation. The contribution from Z — > e + e~ events where both electrons are 
identified as photons, is estimated from simulation. 

The 77 invariant mass distribution is shown in Fig. 5, along with the signal of a 130 GeV Higgs boson. No excess 
over predictions is observed, and limits are set on the possible cross section for the production of such events. The 
limits are also shown in Fig. 5. 

For higher values of rrih, the h — > W + W~ decay mode becomes dominant. The CDF Collaboration has exploited 
this by using a SM Higgs boson search in the channel W h — » W ± W + W~ in a 1.9 fb" 1 dataset, and interpreting its 
result in the context of a search for fermiophobic Higgs bosons [8] . The main difference with the corresponding SM 
search is that the branching fraction for this Higgs boson decay mode is relatively high even for lower values of run- 

The requirement of two isolated leptons of the same charge sign results in a nearly background free sample. The 
remaining backgrounds are due to fake leptons and their rate, after the requirement of isolated tracks with certain 
energy deposits in the hadron calorimeters, is estimated from inclusive jet samples. No signal is observed in the 
distributions of relevant kinematic variables (see e.g. Fig. 6), and also the cross section limits obtained are compared 
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Figure 4: Top: visible mass distribution in the e/j, channel (left) and in the ern + fJ/Th channels (right) in the CDF 
search. Bottom: limits in the (mA, tan/3) plane, for /i < (left) and /i > (right). 
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Figure 5: Left: 77 invariant mass distribution in the DO h — > 77 search. Right: cross section limit as a function of mi,. 



with predictions for SM and fermiophobic Higgs boson production in Fig. 6. 



3. CHARGED HIGGS BOSONS 



Charged Higgs bosons have previously been searched for at the LEP Collider, in their decays to tv or cs final 
states. Assuming that these two decay modes saturate the H^ decays, masses below 79 GeV were excluded. 

If H are relatively light, m(H ± ) < mt — nib, their dominant production mechanism at the Tevatron collider is 
in decays of top quarks, t — v bH + . This decay mode may be an important one either at low or at high values of 
tan/3 (due to the mt cot (3 and mf,tan/3 terms in the coupling, respectively). For moderate to high tan/3 it decays 
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Figure 6: Left: transverse momentum of the dilepton system for the p e T channel in the CDF Wh — > WWW search. 
Right: resulting cross section limits as a function of mh, compared with the predictions for SM and fermiophobic Higgs boson 
production. 

primarily to rP, while for low values the decay is mostly to cs. 

The DO Collaboration has re-interpreted existing ti production cross section measurements in the context of a 
search for the decay H + — > cs under the pessimistic assumption that m(_ff ± ) « mw [9]. The only observable 
consequence of such a scenario is an apparent increased branching fraction for hadronic decays of the W boson; this 
can be tested by a comparison of the measured di-lepton and £+jets tt cross sections. From the observed ratio, 
CT (PP ~~ * ^)i?+iets/ <7 (PP ~ *■ tE)li — l«21_o '26' application of the Feldman-Cousins likelihood ordering principle [10] 
leads to the result B(t -v H+b) < 0.35 at 95% CL. 

A new CDF analysis, on a 2.2 fb^ 1 dataset, has been carried out looking for the same decay but under the 
assumption that m(H ) > M\y [11]. The analysis is restricted to i+jets ti candidate events, requiring in addition 
four jets in the final state, two of which must be b tagged. 

The invariant mass of the remaining two jets is computed. Neglecting backgrounds and resolution effects, the 
resulting distribution should yield the W boson mass. If part of the decays are to charged Higgs bosons, the invariant 
mass distribution is modified. A conrparison between measured and predicted invariant mass distributions is shown 
in Fig. 7, with in addition the expected upper limit at 95% CL for the signal of a putative charged Higgs boson 
of mass 120 GeV. No such deformation of the di-jet invariant mass distribution is observed, and upper limits on 
B(t — » H + b) ■ B(H + — > cs) are also given in Fig. 7. 



4. DOUBLY CHARGED HIGGS BOSONS 



A more exotic extension of the SM Higgs sector is to posit the existence of Higgs multiplets with weak isospin 
I > s, leading to doubly charged Higgs bosons 7J ±± . In general, such extensions are severely constrained by the 
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/m% cos 2 Ow = 1 at tree level if only Higgs doublets are present but can be altered 



p parameter, which is p 

significantly in the presence of I > h multiplets. However, these constraints can be evaded if the doubly charged 
Higgs boson coupling to W bosons vanishes. Two models exist, with "left-handed" ("right-handed") Higgs bosons 
coupling to left-handed (right-handed) fermions, respectively. 

The DO Collaboration has searched in a 1.1 fb -1 dataset for such doubly charged Higgs bosons produced through 
the process pp — ► Z/7* — > H ++ H~~ , with both Higgs bosons decaying to muons, H ±:t — » p^p.^ [12]. As in the case 
of the CDF search for fermiophobic Higgs bosons (see Sect. 2.4), a clean sample is obtained by the requirement of 
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Figure 7: Left: di-jet invariant mass in the CDF t — + H + b, H + — > cs search. The expected deformation due to a 120 
GeV H is also shown, for the purpose of illustration. Right: resulting upper limits on the product of branching fractions 
B(t-> H+b) -B(H+ ->-cs). 



two like-sign isolated muons. Further sensitivity is achieved by requiring a third muon to be present. The dominant 
background is from WZ and ZZ events, with muons associated with a jet passing the isolation criteria. The isolation 
fake rate is known and is applied to scale these backgrounds, inferred from simulation. The number of heavy flavour 
events (with again muons resulting from semi-muonic decays appearing to be isolated) is estimated from the number 
of such events with opposite charge sign. 

Three events remain after all selection criteria, consistent with background expectations as shown in Fig. 8. The 
resulting cross-section limit, obtained as a function of the Higgs boson mass, is compared with the theoretical cross- 
section. This is similarly shown in Fig. 8. Lower limits of 150 GeV(127 GeV) for left-handed (right-handed) doubly 
charged Higgs bosons are obtained. 
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Figure 8: Left: like-sign di-muon invariant mass after all selection criteria in the DO search for doubly charged Higgs bosons. 
The dark- and light-hatched areas represent the expected multi-jet and WZ/ZZ backgrounds; the open histogram represents 
the expected signal due to a left-handed H with a mass of 140 GeV. Right: resulting cross section limits as a function of 
the doubly charged Higgs boson mass. 
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